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Abstract 

We compute the spectrum and several critical amplitudes of the 
two dimensional Ising model in a magnetic field with the transfer 
matrix method. The three lightest masses and their overlaps with the 
spin and the energy operators are computed on lattices of a width 
up to Li = 21. In extracting the continuum results we also take 
into account the corrections to scaling due to irrelevant operators. In 
contrast with previous Monte Carlo simulations our final results are 
in perfect agreement with the predictions of S-matrix and conformal 
field theory. We also obtain the amplitudes of some of the subleading 
corrections, for which no S-matrix prediction has yet been obtained. 
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1 Introduction 



In these last years there has been much progress in the study of 2d spin 
models in the neighbourhood of critical points. The scaling limit of such 
models is described in general by the action A obtained by perturbing the 
conformal field theory (CFT) which describes the critical point with one (or 
more) of the relevant operators which appear in the spectrum of the CFT. 

A = Ao + X J d^x(j){x) (1) 

where is the action of the CFT at the critical point |I[ and (f){x) is the 
perturbing operator. Few years ago A. Zamolodchikov in a seminal paper 
suggested that in some special cases these perturbed theories are equivalent 
to relatively simple quantum field theories whose mass spectrum and S- 
matrix are explicitly known. Later it was realized that these theories had 
a deep connection with the Dynkin diagrams of suitable Lie algebras and, 
from the exact knowledge of the S-matrix, several other informations, and 
in particular some critical amplitudes were obtained (for a review, see for 
instance @]). While these results have formally the status of conjectures, 
they successfully passed in these last years so many tests that they are now 
universally accepted. The most fascinating example of these S-matrix models 
is the Ising model perturbed by an external magnetic field, which is also the 
model which was originally studied by Zamolodchikov in 0. This model is 
highly non-trivial. Its spectrum contains 8 stable scalar particles, all with 
different masses. Both the masses and the entries of the S-matrix are based 
of the numerology of the Eg exceptional Lie algebra. In particular the ratio 
between the first two masses is predicted to be the "golden ratio" 1712/1711 = 
2 cos(|). The simplest realization of this QFT is the 2d Ising model ai (3 = (3c 
in presence of an external magnetic field h. However there are several other 
models which belong to the same universality class. In particular, the first 
numerical check of the predictions of was performed on the Ising quantum 
spin chain in which the first few states of the spectrum were precisely 
observed. Another interesting realization was presented in |^, where the 
dilute ^3 IRF (Interaction Round a Face) model was solved exactly and the 
predicted spectrum of states was found |]^. 

Despite these successes, little progress has been achieved in testing Zamolod- 
chikov's proposal directly in the 2d Ising spin model. Even more, it is exactly 
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for this model that one faces the only existing discrepancy between Zamolod- 
chikov's results and Monte Carlo simulations. 

The Ising spin model in a magnetic field was studied numerically in ^]. 
In both papers the authors studied the spin-spin correlator and did not find 
the spectrum predicted by . Their data were compatible with the presence 
in the spectrum of only the lowest mass state. The explanation suggested 
in 1^, was that probably the higher masses had a negligible overlap am- 
plitude with the spin operator. However, later, in |[T^ these overlaps were 
evaluated explicitly in the S-matrix framework and turned out to be of the 
same order of magnitude as the overlap with the lowest mass state. 

In this paper we shall address this problem. We shall show that Zamolod- 
chikov's proposal (and the calculations of |TD[) is correct also in the case of 



the 2d Ising spin model and that the apparent disagreement was due to the 
fact that it is very difficult to extract a complex spectrum from a multi- 
exponential fit to the spin-spin correlator. We have been prompted to this 
explanation by another example that we recently studied, in which exactly 



the same phenomenon happens: the 3d Ising model [11 1. In this case also 



multi-exponential fit to the spin-spin correlator seems to indicate the pres- 
ence of a single state in the spectrum, while using a suitable variational 
method and diagonalizing a set of improved operators one can clearly see the 
rich spectrum of the model. 

While in previous numerical works p], 3 the model was studied by using 
Monte Carlo simulations we tried in the present paper a different approach 
based on the exact diagonalization of the transfer matrix. 

This approach has various advantages: it gives direct access to the mass 
spectrum of the model and allows to obtain numerical estimates of various 
quantities with impressively small uncertainties. However it has the serious 
drawback that only transfer matrices of limited size can be handled and it 
is difficult to extract from them the continuum limit results in which we are 
interested. During the last years various strategies have been elaborated to 
attack this problem, but all of them are affected by systematic errors whose 
size is in general unknown. 

In this paper we propose a new approach based on the fact that, by using 
the exact solution of the Ising model at the critical point, one can construct 
very precise expansions for the scaling functions in powers of the perturbing 
field. More precisely, thanks to the knowledge of the spectrum of the model, 
it is possible to list all the irrelevant fields which may appear in the effective 
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Hamiltonian and select them on the basis of the symmetry properties of the 
observables under study. 

Our strategy could be summarized as follows. 

• Choose a set of values of h for which the correlation length is much 
smaller than the maximum lattice size that we can studyQ Diagonalize 
the transfer matrix for various values of the transverse size of the lattice 
and extract all the observables of interest. 

• Extrapolate the numbers thus obtained to the thermodynamic limit. 
Thanks to the very small correlation length, the finite size behaviour is 
dominated by a rapidly decreasing exponential and the thermodynamic 
limit can be reached with very small uncertainties (we list in a set of 
tables at the end of the paper the results that we obtained in this way). 

• Construct for each observable the scaling function keeping the first 7 
or 8 terms in the expansion in powers of the perturbing field. 

• Fit the data with these truncated scaling functions. By varying the 
the number of input data and of subleading terms used in the scaling 
functions we may then obtain a reliable estimate of the systematic 
deviations involved in our estimates (see the discussion in sect. 6). 

Which are the observables of interest mentioned above? 

Usually, when looking at the scaling regime of statistical models one 
can study only adimensional amplitude ratios which are the only quanti- 
ties which, thanks to universality, do not depend on the details of the lattice 
models, but only on the features of the underlying QFT. However the Ising 
model can be solved exactly at the critical point also on the lattice and ex- 
plicit expressions for the spin-spin and energy-energy correlators are known. 
This allows to write the exphcit expression in lattice units of the amplitudes 
evaluated in the framework of the S-matrix theory. Thus one is able to pre- 
dict not only adimensional amplitude ratios but also the values of the critical 
amplitudes themselves. This greatly enhances the predictive power of the S- 
matrix theory and makes much more stringent the numerical test that we 
perform. 

^In particular we decided to keep the ratio < 0.1. This means that we only studied 
values of h for which the correlation length was smaller than two lattice spacings. 
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The final result of our analysis is that all the observables that we can 
measure perfectly agree with the S-matrix predictions. 

In particular we obtain very precise estimates for the first three masses, 
for several critical amplitudes and, what is more important, for the overlap 
amplitude of the first two masses with the spin and energy operators, a result 
which had never been obtained before. 

We also measure the amplitude of some of the subleading corrections in 
the scaling functions, for which no S-matrix prediction exists for the moment. 
In particular we found that the amplitude of the corrections due to the energy 
momentum tensor in translationally invariant observables is compatible with 
zero. 

This paper is organized as follows. In sect. 2 we introduce the model in 
which we are interested, collect some known results from S-matrix theory 
and finally give the translation in lattice units of the critical amplitudes 
evaluated in the S-matrix framework. In sect. 3 we construct the scaling 
functions. This only requires the use of very simple and well known results 
of Conformal Field Theory. Notwithstanding this it turns out to be a rather 
non trivial exercise. Since it could be a result of general utility (it could be 
extended to other models for which the CFT solution is known or to other 
quantities of the Ising model in a magnetic field that we have not studied in 
the present paper) instead of simply giving the results, we derived the scaling 
functions explicitly and tried to give as much details as possible. Sect. 4 is 
devoted to a description of the transfer matrix method. In sect. 5 we deal 
with the thermodynamic limit while in sect. 6 we analyze the transfer matrix 
results and give our best estimates for the critical amplitudes in which we 
are interested. Sect. 7 is devoted to some concluding remarks. 

To help the reader to reproduce our analysis (or to follow some alternative 
fitting procedure) we list in four tables at the end of the paper the data that 
we obtained with the transfer matrix approach. 

2 Ising model in a magnetic field 

In this section we shall review the existing theoretical informations on the 
Ising model in a magnetic field. This will require four steps. First (in 
sect. 2.1) we shall define the lattice version of the model, discuss its action 
and define the observables in which we shall be interested in the following. 
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Then (in sect. 2.2) we shall turn to the continuum version of the theory, 
described in the present case by the action: 

Ao + h j d^x (t{x) (2) 

where (t{x) is the perturbing operator. In particular we shall discuss, within 
the framework of the renormalization group, the expected scaling behaviour 
of the various quantities of interest and define the corresponding critical 
amplitudes. In sect. 2.3 we shall use the knowledge of the S- matrix of the 
model to obtain the value of some of the amplitudes of interest by using 
the Thermodynamic Bethe Ansatz (TBA) and the form factor approach. 
Finally in sect. 2.4 wc shall turn back to the lattice model and show how the 
continuum results can be translated in lattice units. 



2.1 The lattice model 

The lattice version of the Ising model in a magnetic field is defined by the 
partition function 

where the field variable cr„ takes the values {±1}; n = {uq, ni) labels the sites 
of a square lattice of size Lq and Li in the two directions and (n, m) denotes 
nearest neighbour sites on the lattice. In our calculations with the transfer 
matrix method we shall treat asymmetrically the two directions. We shall 
denote Uq as the "time" coordinate and Ui as the space one. The number of 
sites of the lattice will be denoted by = LqLi. In the thermodynamic limit 
both Lq and Li must go to infinity and only in this limit we may recover the 
results of the continuum theory. In our actual calculations with the transfer 
matrix method we shall study finite values of Li and then extrapolate the 
results to infinity. This extrapolation induces systematic errors which are 
the main source of uncertainty of our results, since the rounding errors in the 
transfer matrix diagonalization are essentially negligible. In sect. 6 below, 
we shall discuss these systematic errors and estimate their magnitude. 

In order to select only the magnetic perturbation, the coupling /3 must 
be fixed to its critical value 

= = 1 log (\/2 + 1) = 0.4406868... 
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by defining hi = (3cH we end up with 

Z{hi)= ^ e^"^<"^'">'^™'^"'^'''^"'^" . (4) 

hi denotes the lattice discretization of the magnetic field h which appears 
in the continuum action eq. (0). It must be, for symmetry reasons, an odd 
function of h. 

2.1.1 Lattice operators 

It is useful to define the lattice analogous of the spin and energy operators 
of the continuum theory. They will correspond to linear combinations of 
the relevant and irrelevant operators of the continuum theory with suitable 
symmetry properties with respect to the Z2 symmetry of the model (odd for 
the spin operator and even for the energy one). Near the critical point this 
linear combination will be dominated by the relevant operator and the only 
remaining freedom will be a conversion constant relating the continuum and 
lattice versions of the two operators (we shall find this constants in sect. 2.5). 
The simplest choices for these lattice analogous are 

• Spin operator 

ai{x) = (T^ (5) 

i.e. the operator which associates to each site of the lattice the value 
of the spin at that site. 

• Energy operator 




where the sum runs over the four nearest neighbour sites y of x. 
represents a constant "bulk" term which we shall discuss below. 

The index I indicates that these are the lattice discretizations of the contin- 
uous operators. We shall denote in the following the normalized sum over all 
the sites of these operators simply as 
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2 Observables 

► Free Energy 

The free energy is defined as 

f{hi)^j^\ogiZihi)) . (8) 

It is important to stress that f{hi) is composed by a "bulk" term fb{hi) 
which is an analytic even function of hi and by a "singular" part fs{hi) 
which contains the relevant informations on the theory as the critical 
point is approached. The continuum theory can give informations only 
on fg. The value of /^(O) can be obtained from the exact solution of 
the lattice model at hi = 0, P = Pc (see [p!2| ) 



2G 1 

fb = — + -log 2 = 0.9296953982... (9) 

TT 2 

where G is the Catalan constant. 
Magnetization 

The magnetization per site M{hi) is defined as 

M{hi)^j^^{\og Z{hi)) = j^{Y:cr,) . (10) 
Hence we have 

M{hi) = {ai) . (11) 

Magnetic Susceptibility 

The magnetic susceptibility x is defined as 

Internal Energy 

We define the internal energy density E{hi) as 

^(^0 = ^(E • (13) 

(n,m) 
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As for the free energy, also in this case one has a bulk analytic contri- 
bution Eb{hi) which is an even function of hi. Let us define = Eb{0). 
The value of -^^(O) can be easily evaluated (for instance by using 
Kramers- Wannier duality) to be = Let us define E{hi) = 

E{hi) — eb, we have 

E{hl) = ^( E ^n^m) - \ . (14) 
{n,m) ^ 



Hence we have 



Eihi) = {ei) . (15) 



As usual the internal energy can also be obtained by deriving the free 
energy with respect to /3. However it is important to stress that, due 
to the magnetic perturbation (see eq.(^) in performing the derivative 
we also extract from the Boltzmann factor a term proportional to Hai. 
Hence we have: 

This observation will play an important role in the following. 
2.1.3 Correlators 

We are interested in the spin-spin and in the energy-energy connected corre- 
lators defined as 

G.,,(r) = {cJiiO)aiir)) - {a^f = (aKO)aKr)), , (17) 

G,,(r) = (Q(0)Q(r))-(Q)^ = (Q(0)eKr)). . (18) 

For a nonzero magnetic field these correlators are very complicated, unknown, 
functions of h and r, however a good approximation in the large distance 
regime r ^ oo isQ 



TT 



^For a discussion of the limits of this approximation and of the corrections which must 



be taken into account when the short distance regime is approached see [13 
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where the sum is over the low laying single particle states of the spectrum, 
mi{h) denotes their mass the functions F^{h) their overlap with the a oper- 
ator. Similarly we have 

2^ . (20) 

where the spectrum is the same as for the spin-spin correlator but the overlap 
constants are different. 

A particular role is played by the lowest mass nii which gives the dom- 
inant contribution in the large distance regime. Its inverse corresponds to 
the (exponential) correlation length ^ of the model and sets the scale for 
all dimensional quantities in the model. In particular the "large distance 
regime" mentioned few lines above means "large with respect to the correla- 
tion length" . 

In the approximation of eqs. (|19|) and (^0|) one is neglecting the cut-type 
contributions which appear above the two-particle threshold i.e. at twice the 
value of mi. For this reason we shall concentrate in the following only on 
the three first states of the spectrum which are the only ones which lie below 
such threshold (see eq. ( ^D] ) below). 

2.1.4 Time slice correlators 

It is very useful to study the zero momentum projections of the above de- 
fined correlators. They are commonly named time slice correlators. The 
magnetization of a time slice is given by 

^no = X! '^{no,ni) ■ (21) 

The time slice correlation function is then defined as 

GlAr) = E{(^no^no+.) - {SrS) (22) 

no 

where the index indicates that this is the zero momentum projection of the 
original correlator. Starting from eq. ([T9|) it is easy to show that in the large 
r limit G°^(r) behaves as 

'^LW _ ^ ]^7Wl^g-m,(h)|r| _ ^23) 
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A similar result, with the obvious modifications, holds also for G° 



2.2 Critical behaviour 

In this section we discuss the critical behaviour of the model by using stan- 
dard renormalization group methods, keeping in the expansions only the first 
order in the perturbing field. Both the results and the analysis are well known 
and can be found in any textbook. We report it here since it will serve us 
as a starting point for the more refined analysis which we shall perform in 
sect. 3 below. 



2.2.1 Critical indices 



The starting point of the renormalization group analysis is the singular part 
of the free energy fs(t,h) (where t is the reduced temperature). Standard 



renormalization group arguments (see for instance [jT4|) allow to write fs in 
terms of a suitable scaling function $: 



fs{t,h) 



Uh 



d/Vh 



[\uh/uh,\y^/y\ 



(24) 



where Utg and are reference scales that depend on the model. Uh and Ut 
denote the scaling variables associated to the magnetic and energy operators 
respectively and i/h, yt are their RG-exponents. Ut and Uh do not exactly 
coincide with t and h but are instead analytic functions of them. The only 
constraint is that they must respect the Z2 parity of t and h . Near the 
critical point we may suitably rescale $ so as to identify Ut = t and Uh = h. 
Thus, setting t = we immediately obtain the asymptotic critical behaviour 
of/. 

fsOc\h\''/y'^ . (25) 

Taking the derivative with respect to h (or t) and then setting t = we 
can obtain from eq. ( ^^ also the asymptotic critical behaviour of the other 
observables in which we are interested 



X oc \hfy^-^ 



(26) 
(27) 
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E oc I /i I ('^-f . (28) 

From the exact solution of the Ising model at the critical point we know 
that Uh = Y Ut = ^- Inserting these values in the above expressions we 
find 

/s(x|/z|i (29) 
M(x\h\T^ (30) 
Xoc|/i|-^ (31) 
Eoc\h\TB . (32) 
The masses rrii have as scaling exponent, as usual, l/yn, hence 

TTij oc . (33) 

Finally from the definitions of eqs. (|19D and (|2^) we see that the overlap 
amplitudes behave as adimensional constants. 



2.2.2 Critical amplitudes 

In order to describe the scaling behaviour of the model we also need to know 
the proportionality constants in the above scaling functions. These constants 
are usually called critical amplitudes. Using the results collected in eqs. (pQ])- 
( |32D we have the following definitions: 



Af = lim / /i 15 , Am = lim M h , Ay = lim x ^^'^ , (34) 

h—*0 ?i— >0 h^O 
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Ae = limE h w , A^- = lim mj h w , (35) 

Ap. = lim , Af^ = lim F' . (36) 

Notice for completeness that in the literature (see for instance [0) the 
amplitudes A^, Am and A^^ are usually denoted as 

= F, , Am = Dc^^ , A^, = y . (37) 
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We shall show in the next section that all these amplitudes can be exactly 
evaluated in the framework of the S'-matrix approach. As a preliminary step 
let us notice that since M and x ci-re obtained as derivatives of / we have 

Am = —Af , A^ = —Am • (38) 



2.2.3 Universal amplitude ratios 

From the above critical amplitudes one can construct universal combinations 
which do not depend on the particular realization of the model. For this 
reason they have been widely studied in the literature. In particular there are 
two "classical" amplitude combinations which involve the critical amplitudes 



defined above (see for instance ||T5[). They are: 

= VD,B'^ Q2 = (F/Fe) iU^o) ^ (39) 

where we used the notations of eq. (^). F and denote the critical ampli- 
tudes of the susceptibility and exponential correlation length for h = and a 
small positive reduced temperature, while B denotes the critical amplitude 
of the magnetization for h = and a small negative reduced temperature. 
Notice however that, since (as we mentioned in the introduction) we are 
able to give the explicit relation between lattice and continuum expectation 
values, we are not constrained to study only universal combination but can 
determine exactly the various critical amplitudes. 



2.3 S-matrix results 

In 1989 A. Zamolodchikov suggested that the scaling limit of the Ising 
Model in a magnetic field could be described by a a scattering theory which 
contains eight different species of self-conjugated particles Aa, a = 1, . . . , 8 
with masses 



m2 


= 2mi cos 


TT 

5=' 


1.6180339887..) mi. 




= 2mi cos 


TT 

30 ~ 


(1.9890437907..) mi. 


m4 


= 2m2 cos 


7n 

30 ~ 


(2.4048671724..) mi. 
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ms = 2m2 cos ^ = (2.9562952015..) mi, (40) 

me = 2m2 cos ^ = (3.2183404585..) mi, 

n 7n 

my = 4m2 COS — COS — = (3.8911568233..) mi , 
5 30 

7r 27r 

mg = 4m2 cos — cos — = (4.7833861168..) mi 
5 15 

where mi{h) is the lowest mass of the theory. As mentioned above it coincides 
with the inverse of the (exponential) correlation length. Few years later, from 
the knowledge of the S-matrix of the theory V. Fateev [|l6l obtained explicit 
predictions for some of the critical amplitudes defined above. 

In order to evaluate the amplitudes one must first fix the normalization of 
the operators involved which can be set, for instance, by fixing the constant 
in front of the long distance behaviour of the correlators at the critical point. 
It is important to make explicit this normalization choice, since it will allow 
us, by comparing with the corresponding correlators in the lattice theory to 
convert explicitly the continuum results in lattice units. Following the choice 
of [jl6[ we assume: 



\x 



\x\ 



oo 



(41) 



(e(x)e(O)) 



\x\ 



With these conventions one finds 16 



where 



A 



mi 



c 



A, 



\x\ 



oo. 



(sin^ + sin^ + sin^) 



(42) 

(43) 
(44) 



C 



4 sin IF (I) /4vr^r(|)F^(l|)y 



r(i)r(^) V r(i)n(A) 

From Af one immediately obtains Am and A^^. 



4.40490858... 



(45) 
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The amplitude Ae requires a more complicated analysis. Its exact ex- 
pression has been obtained only recently in [|17| 

Ae = 2.00314... . (46) 

We summarize in tab. |l| these S-matrix predictions for the critical ampli- 
tudes. 



A 


= 4.40490858.. 




= 1.19773338.. 


Am 


= 1.27758227.. 


\ 


= 0.08517215.. 


Ae 


= 2.00314.. 



Table 1: Critical amplitudes. 

From these critical amplitudes, and using the values of F, B and one 
immediately obtains the classical amplitude ratios defined above (see for 
instance |jl8[). They are reported in tab. ^ . 



= 6.77828502.. 
Qa = 3.23513834.. 



Table 2: Classical amplitude ratios. 

Finally, the critical overlap amplitudes Ae<^ and Ae^ were evaluated in 
[p!0| , p!9[] . They are reported in tab. |^ and § . 

2.4 Conversion to lattice units 

While the values listed in tab. |^, ^ and ^ are universal, the amplitudes listed 
in tab. |I| depend on the details of the regularization scheme. Thus some 
further work is needed to obtain their value on the lattice. We shall denote 
in the following the lattice critical amplitudes with an index I. Thus, for 
instance, 

A'm = lim>^) h;^^ , (47) 

/l7->0 
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Ap,r = 
^1 


-0.64090211.. 


Ap" = 

^2 


0.33867436.. 




-0.18662854.. 


Ap" = 


0.14277176.. 


Ap<r = 


0.06032607.. 


Apcr = 


-0.04338937.. 


Ap. = 


0.01642569.. 


Ap-r = 


-0.00303607.. 



Table 3: Critical overlap amplitudes for the spin operator. 



Ape = 


-3.70658437.. 


Ape = 


3.42228876.. 


Ape = 


-2.38433446.. 


Ape = 


2.26840624.. 


Ape = 


1.21338371.. 


Ap: = 

Ape = 


-0.96176431.. 


0.45230320.. 


Ape = 


-0.10584899.. 



Table 4: Critical overlap amplitudes for the energy operator. 

to be compared with the continuum critical amplitude defined in eq. (p^ 

Am = lim((7) h'^ . (48) 

In order to relate the lattice results with the continuum ones we must 
study the relationship between the lattice operators and the continuum ones. 
In general the lattice operators will be given by the most general combina- 
tion of continuum operators compatible with the symmetries of the lattice 
operator multiplied by the most general analytic functions of t and h (with a 
parity which is again constrained by the symmetry of the operators involved). 
Thus, for instance, anticipating the discussion that we shall make in sect. 3, 
we have 

ai = f^it,h)a + f,{t,h)<Pi (49) 
where fo{t, h) and /j(t, h) are suitable functions of t and h and with 0j we 
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denote all the other fields of the theory (both relevant and irrelevant) which 
respect the symmetries of the lattice. 

A similar relation also holds for the energy operator: 

ei = g'o{t,h)e + g,{t,h)^i . (50) 

Finally, also hi is related to the continuum magnetic field h hj a relation 
of the type 

hi = boit,h)h (51) 

where bolt, h) must be an even function of h . 

At the first order in t and h these combinations greatly simplify and es- 
sentially reduce to a different choice of normalization between the continuum 
operators and their lattice analogous: 

ai = R^a , ei = R^e , hi = Rhh (52) 

where R^, Re and Rh are three constants which correspond to the h ^ 0, 
t — > limit of the /q , and bo functions. 

If we want to compare the S-matrix results discussed in the previous sec- 
tion with our lattice results we must fix these normalizations^. The simplest 
way to do this is to look at the analogous of eqs. (P1^2|) at the critical point 
(namely for hi = 0) [0 . 

In fact, if hi = it is possible to obtain an explicit expression for the 
spin-spin and energy-energy correlators (for a comprehensive review see for 
instance [^) directly on the lattice, for any value of (3. Choosing in par- 



ticular (3 = (3c, and looking at the large distance behaviour of these lattice 
correlators we may immediately fix the normalization constants. Let us look 
first at Ra- 

We know from ||2^ that: 

R^ 

{(riCXj)h=o = — ^ (53) 

where rij denotes the distance on the lattice between the sites i and j and 

^2 _ g3r(-i)25/24 _ 7Q33g^^^ ^ (^54) 

■^This essentially amounts to measure all the quantities in units of the lattice spacing. 
For this reason we can fix in the following the lattice spacing to 1 and neglect it. 
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By comparing this result with eq. (0) we find 



= 0.83868... 



(55) 



From this we can also obtain the normalization of the lattice magnetic 
field which must exactly compensate that of the spin operator in the pertur- 
bation term ha. We find: 



Rh 



(Ra 



,-1 



1.1923... 



(56) 



Combining these two results we obtain the value in lattice units of the 
constant A^^ 

A\, = {R^Y^'^^'Am = 1.058... . (57) 

From this one can easily obtain also Aj, and Am^. 

Let us look now at R^. In the case of the energy operator the connected 
correlator on the lattice, at /i; = and for any value of {3 has the following 
expression 



23 



(Q(O)Q(r)). 



'5^ 



TT 



Kl{5r) - Kl{5r) 



(58) 



where Kq and Ki are modified Bessel functions, 5 is a parameter related to 
the reduced temperature, defined as 



5 = A\(3-(3, 



(59) 



and with the index c we denote the connected correlator (notice that thanks 
to the definition (|T^ no disconnected part must be subtracted at the critical 
point and the index c becomes redundant). This expression has a finite value 
in the 5^0 limit (namely at the critical point). In fact the Bessel functions 
difference can be expanded in the small argument limit as 



Kl{5r) ~ Kl{5r) 



thus giving, exactly at the critical point: 



1 



(5r) 



+ ... 



(60) 



(Q(O)Q(r)) 



7rr 



(61) 
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By comparing this result with eq. ( ^21) we find 



Re = - 
TT 

and from this we obtain the expression in lattice units of 
A[ = {R„f'^\R,)A, = 0.58051... . 
Our results are summarized in tab. |^ . 



A 

A'- 

A'. 



mi 



M 



4.01039911... 
0.99279949... 
1.05898612... 
0.07059907... 
0.58051... 



(62) 



(63) 



Table 5: Critical amplitudes in lattice units. 



2.4.1 Alternative derivation of R^ 

In this section we discuss, for completeness, an alternative derivation of R^. 
It can be used in those cases in which the correlators are not known, but the 
internal energy is known on a finite size lattice at the critical point. Then 
i?£ can be obtained by comparing the finite size behaviour of the internal 
energy on the lattice with that predicted by conformal field theory in the 
continuum. In the case of the Ising model, thanks to the beautiful work by 
Ferdinand and Fisher |]12|, we know that on a square lattice of size Lq x Li 



with Lq > Li with periodic boundary conditions the internal energy must 
scale as: 

where 'i9i(r) denotes the i*^ Jacobi theta function and t = ij^. 

The same behaviour can be studied in the continuum theory, by using 
CFT techniques. The result is 



A{t) + Mr) + Mr) L, ■ ^ ^ 
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By using the relation 



Mr)' = vr Mr)Mr)Mr) (66) 

which allows to express the derivative of the "t^i (r) in terms of ordinary theta 
functions we see that the two equations (|64D and (|65D agree only if we choose, 
as we did in the previous section, R^ = K 

3 Scaling functions 

In this section we shall construct the scaling functions for the various quan- 
tities in which we are interested. Our aim is to give the form (i.e. the value 
of the scaling exponents) of the first 7-8 terms of the expansion in powers of 
h of the scaling functions and at the same time to identify the operators in 
the lattice Hamiltonian from which they originate. To this end we shall first 
deal in sect. 3.1 with the theory at the critical point. We shall in particular 
discuss its spectrum, which can be constructed explicitly by using CFT tech- 
niques. Next, in sect. 3.2, we discuss in the framework of the renormalization 
group approach the origin of the subleading terms in the scaling functions, 
and show how to obtain their exponents from the knowledge of the renor- 
malization group eigenvalues Ui of the irrelevant operators. While in general 
this analysis is only of limited interest since the Ui of the irrelevant operators 
are unknown, in the present case, thanks to the CFT solution discussed in 
sect. 3.1, it becomes highly predictive and will allow us to explicitly construct 
in sects. 3.3 and 3.4 the scaling functions. In particular in sect. 3.3 we shall 
list all the irrelevant operators which may appear in the effective Hamiltonian 
and discuss their symmetry properties, while in sect. 3.4 we shall write the 
scaling functions and identify the operators involved in the various scaling 
terms. 

3.1 The Ising model at the critical point 

The Ising model at the critical point is described by the unitary minimal 
model with central charge c = 1/2 [|l|. Its spectrum can be divided into 
three conformal families characterized by different transformation properties 
under the dual and Z2 symmetries of the model. They are the identity, spin 
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and energy families and are commonly denoted as [/] , [a] , [e] . Let us discuss 
their features in detail. 

• Primary fields 

Each family contains a relevant operator which is called primary field 
(and gives the name to the entire family). Their conformal weights 
are hj = 0, = 1/16 and = 1/2 respectively. The relationship 
between conformal weights and renormalization group eigenvalues is: 
y — 2 — 2h. Hence the relevant operators must have h < 1. 

• Secondary fields 

All the remaining operators of the three families (which are called sec- 
ondary fields) are generated from the primary ones by applying the 
generators L_j and L_j of the Virasoro algebra. In the following we 
shall denote the most general irrelevant field in the [a] family (which 
are odd with respect to the Z2 symmetry) with the notation cij and 
the most general fields belonging to the energy [e] or to the identity 
[I] families (which are Z2 even) with e, and rji respectively. It can be 
shown that by applying a generator of index k: L^k or L^k to a field 
(f) (where (p =, /, e, a depending on the case) of conformal weight we 
obtain a new operator of weight h = + k. In general any combina- 
tion of L_j and L_j generators is allowed, and the conformal weight of 
the resulting operator will be shifted by the sum of the indices of the 
generators used to create it. If we denote by n the sum of the indices 
of the generators of type L_j and with n the sum of those of type L_i 
the conformal weight of the resulting operator will be + n + n. The 
corresponding RG eigenvalue will he y = 2 — 2h^ — n — n, hence all the 
secondary fields are irrelevant operators. 

• Nonzero spin operators 

The secondary fields may have a non zero spin, which is given by the 
difference n — n. In general one is interested in scalar quantities and 
hence in the subset of those irrelevant fields which have n — n. However 
on a square lattice the rotational group is broken to the finite subgroup 
C4 (cyclic group of order four). Accordingly, only spin 0,1,2,3 are 
allowed on the lattice. If an operator of the continuum theory has 
spin j G N, then its lattice discretization (pi behaves as a spin j (mod 4) 
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operator with respect to the C4 subgroup. As a consequence all the 
operators which in the continuum limit have spin j = AN with N non- 
negative integer can appear in the lattice discretization of a scalar field. 
This will play a major role in the following. 

• Null vectors 

Some of the secondary fields disappear from the spectrum due to the 
null vector conditions. This happens in particular for one of the two 
states at level 2 in the a and e families and for the unique state at 
level 1 in the identity family. From each null state one can generate, 
by applying the Virasoro operators a whole family of null states hence 
at level 2 in the identity family there is only one surviving secondary 
field, which can be identified with the stress energy tensor. 

• Secondary fields generated by L_i 

Among all the secondary fields a particular role is played by those gen- 
erated by the L_i Virasoro generator. L_i is the generator of trans- 
lations on the lattice and as a consequence it has zero eigenvalue on 
translational invariant observables. Another way to state this results is 
to notice that L_.i can be represented as a total derivative, and as such 
it gives zero if applied to an operator which can be obtained as the in- 
tegral over the whole lattice of a suitable density (i.e. a translationally 
invariant operator). 



3.2 RG analysis for /i 7^ 

We shall discuss the higher order corrections to the RG analysis of sect. 2.2 



along the lines of p5|, to which we refer for a more detailed discussion. The 



only improvement that we make with respect to 1 25] is in the part devoted 
to the contribution due to the irrelevant operators, in which we shall make 
use of the results discussed in the previous section. 

We expect three types of corrections to the asymptotic results reported 
in sect. 2.2: 

a] Analytic corrections. 

They are due to the fact, already mentioned in sect. 2.2, that the actual 
scaling variables in the RG approach are not hi and t but Uh and Ut 
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which are in principle the most general analytic functions of hi and 
t which respect the Z2 parity of hi and t. Let us write the Taylor 
expansion for Uh and Ut, keeping only those first few orders that are 



needed for our analysis (we use the notations of 

Uh = hi[l + Cht + dht" + ehh^i + 0(^^^/^f)] (67) 

u^ = t + hthi + ctt^ + dtt^ + Ctthf + /^/if + 0{t'^,t'^hf) (68) 

The corrections induced by the higher terms in Uh and Ut are of three 
types. 

• The first one is very simple to understand. It is due to the higher 
powers of hi contained in Uh which lead to corrections to the power 
behaviours listed in sect. 2.2.1 which are shifted by even integer 
powers of hi. For instance in the free energy, as a consequence of 
the ehhf term in Uh, we expect a correction of this type: 

/,(/iO = 4i/i/i^(i + 4,3i/id' + ••••) ■ (69) 

with ylj3 = ^Ch- The indices /, 3 in Af^^ only denote the fact 
(that we shall discuss in detail in the next section) that this term 
is the third term in the hi expansion of the scaling function of the 
singular part of the free energy. 

• The second type of correction is due to the terms that depend on 
hi which appear in Ut- Their peculiar feature is that, even if they 
are originated by analytic terms in the scaling variables, they lead 
in general to non analytic contributions in the scaling functions 
For instance, as a consequence of the bthf term in Ut, we find in 
the free energy a correction of the type: 

/,(/iO = 4i/i#(i + 4,2i/^/l'"^ + ••••) • (70) 

with A,, = %^b, and 2 - 1^ = §. 

• The corrections of the third type only appear when studying the 
internal energy. They are due to the terms linear in t which are 
present in Uh and Uf. The most important of these contributions 
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is the one due to the Cht term in Uh which gives a correction pro- 

_8_ 

portional to h^^ to the dominant scahng behaviour of the internal 
energy. We shall discuss these terms in sect. 3.4.3 below. 

b] Corrections due to irrelevant operators in the lattice Hamilto- 
nian. 

These can be treated within the framework of the RG as follows. Let 
us study as an example the case of an irrelevant operator belonging to 
the Identity family. Let us call the corresponding scaling variable M3 
and its RG eigenvalue (since M3 is irrelevant, < 0). In this case 
the dependence of U3 on t and hi is 

= ul + at + bhf -\ . (71) 

Let us for the moment neglect higher order terms and assume ^3 = ^3. 
Then looking again at the singular part of the free energy we find 

/.(t,/.0 = IV/^or/^'^$(^^,«^|/^/|l^^l/^'0 • (72) 

Since Mgl/i^l'^^'/^'' is small as hi ^ it is reasonable to assume that we 
can expand fs in a Taylor series of Ugl/i;!'^^'/^'' (notice that in eq. (|7^ ) fs 
is not singular since it is evaluated at \hi\ > 0). Hence we find (setting 
again t = 0) 

fs = Ihif/^'ia, + a2M°|/ii|l^^l/^'^ + ■ ■ ■) (73) 
where ai, 02, M3 are non- universal constants. 

This analysis can be repeated without changes for any new irrelevant 
operator: M4, 7/4 and so on. As a last remark, notice that on top of these 
non analytic corrections we also expect analytic contributions due to 
the higher order terms contained in eq. ([7I|). 

While in general this analysis is only of limited interest since the yi of 
the irrelevant operators are usually unknown, in the present case we 
may identify the irrelevant operators with the secondary fields discussed 
in 3.1 and use the corresponding RG-exponents as input of our analysis. 

^In general for the irrelevant operators there is no need to tune U3 to to approach 
the critical point. However we shall see below that, for symmetry reasons, u§ = for all 
the irrelevant operators belonging to the [a] and [e] families. 
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c] Logarithmic corrections. 

As it is well known, the specific heat of the 2d Ising model at hi = ap- 
proaches the critical point with a logarithmic singularity. This means 
that in the free energy there must be a term of the type ^ouflog{ut). 
While in general we could expect $o to be a generic function of the ratio 
Ut/v^fl^^ 1 the absence of leading log corrections in M and x strongly con- 
straints this function which is usually assumed to be a simple constant. 
Notwithstanding this, the presence of terms that depend on hi in ut 
implies that log type contributions may appear also in the case t = 0, 
/i; 7^ in which we are interested. These can be easily obtained by 
inserting eq.(|68|) into ^ou^log{ut) and then making the suitable deriva- 
tives and limits [^. In the case of the free energy one obtains a term 
proportional to hflog{hi) which is too high to be observed in our fits. 
However for the internal energy the first contribution is proportional 
to a smaller power of hf. hflog{hi) and must be taken into account in 
the scaling function. 

3.3 The effective lattice Hamiltonian 

Let us call He ft the Hamiltonian which describes the continuum theory at 
the critical point. The perturbed Hamiltonian in the continuum is given 
by: 

H = HcFT + ha . (74) 

The aim of this section is to construct the lattice analogous (which we shall 
call Hiat) of H. 

Notice that Hiat is different from the microscopic Hamiltonian which ap- 
pears in the exponent of eq. (H). Eq. (^ describes the model at the level of 
the lattice spacing. We are instead interested in the large distance effective 
Hamiltonian which one obtains when the short range degrees of freedom are 
integrated out, i.e. after a large enough number of iterations of the Renormal- 
ization Group transformation has been performed. Hiat will contain all the 
irrelevant operators which are compatible with the symmetries of the lattice 
model. In this section we shall first discuss the relation between the lattice 

^Here we follow the convention usually adopted in conformal field theory. In the stan- 
dard notation of classical statistical mechanics one would denote this quantity "Hamilto- 
nian density" rather than "Hamiltonian" . 
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and the continuum operators, then we shall construct the lattice Hamiltonian 
in the hi = case and finally we shall extend our results to the hi ^ case. 

3.3.1 Relations between lattice and continuum operators. 

The lattice operators are given by the most general combination of continuum 
operators compatible with the symmetries of the lattice operator multiplied 
by the most general analytic functions of t and hi (with a parity which is again 
constrained by the symmetry of the operators involved). In the following, to 
avoid a too heavy notation, we shall neglect the t dependence^. 
For the spin operator we have 

c^i = fo{hi)a + hif^{hi)e + f:{hi)a, + hin{hi)e., + hifl{hi)r^,, i e N (75) 

where fi{hi) flihi) and flihi) are even functions of hi. 
For the energy operator we have 

ei = gl{hi)e + hig^{hi)cr + hig^{hi)a, + gl{hi)ti + h^g'{hi)Vi, « ^ N (76) 

where again g^hi) glihi) and gl{hi) are even functions of hi and the h^ term 
in front of gl {hi) is due to the change of sign of the e operator under duality 
transformation at /i/ = (see the discussion at the beginning of sect. 3.3.2). 

Among all the possible irrelevant fields only those which respect the lattice 
symmetries (i.e. those of spin {mod 4)) are allowed in the sums. At this 
stage also irrelevant operators containing L_i or L„i appear in the sums. It 
is only when these operators are applied on translationally invariant states 
(i.e. on the vacuum) that they disappear. This will happen for instance 
when we shall study the mean value of the free energy. 

3.3.2 Construction of Hiat{hi = 0) 

In this case all the operators belonging to the [a] family are excluded due to 
the Z2 symmetry. Also the operators belonging to the [e] family are excluded 
for a more subtle reason. The Ising model (both on the lattice and in the 
continuum) is invariant under duality transformations while the operators 

^The t dependence in the scaling variables of the irrelevant operators plays a role only 
in the construction of the scaling function for the internal energy and we shall resume it 
in sect. 3.4.3 below. 
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belonging to the [e] family change sign under duality, thus they also cannot 
appear in Hiat{hi = 0). Thus we expect 



Hlat = He FT + U^iVi, Vi^ [I] , (77) 

where the are constants. There are however further restrictions: 

• Hlat is a scalar density, hence only operators (p with angular momentum 
j = 4fc, = 0, 1, 2 ■ ■ ■ are allowed. 

• The operator which acts on the Hilbert space of the theory is the space 
integral of the Hamiltonian (density) Hiat- As such it is translational 
invariant and only operators (j) which do not contain the generators L_i 
or L_i survive the integration of eq. (7?) over the space. 



Let us list in order of increasing conformal weight the first few operators 
which fulfill all the constraints: 

<Po = I, 0i = L_2l-2/, 02 = (^-2)'/, h = {L-2?I, (78) 

04 = L_4l, 05 = Z-4/, 06 = 1^31^3! ... . (79) 

Some of these fields have a natural interpretation. 0o gives rise to the bulk 
contribution in the free energy (see the discussion of sect. 2.1.2). 0i, 02, 03 
are related to the energy momentum tensor: TT, T^,T^ respectively. All 
the fields listed above except the identity and 06 have the same conformal 
weight h(f, = 4. The corresponding RG eigenvalue is = —2. The field 06 
has conformal weight h^^ = 6 and RG eigenvalue y^^ = —4. 



3.3.3 Extension to hi ^ 

Mimicking the continuum case we have, also on the lattice, 

Hiatihi) = Hiatihi = 0) + hiai . (80) 
Inserting the expression of ai of eq. (|7^) we find 

Hiatihi) = HcFT + Ui{hi) (81) 

where this time there is no more restriction coming from the Z2 symmetry 
and duality, hence 0j denotes here the most general operator of the spectrum 
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with spin j = 4k, k = 0,1,.. The Ui{hi) are even or odd functions of hi, 
depending on the parity of (pi but in the even sector only for the operators 
belonging to the identity family Ymiin^QUiihi) ^ (according to eq. ( [7T| ) we 



have Ui{hi = 0) = m° ). For the operators belonging to the energy family the 
first nonzero contribution in the Ui{hi) functions is of order hf. 

Let us list, starting from those with the lowest conformal weight, the new 
operators which were not present in Hiat{hi = 0). For future convenience 
let us separate those which do not contain generators from the 

remaining ones. 

A] Operators which are not generated by 

• Operators belonging to [a] 

In the [a] family the lowest ones are L^^a, L^^a and L^^L^a. In 
fact L_i(T disappears for translational invariance and due to the 
null vector equation the L„20" operator which appears at level 2 
can always be rewritten as L'^^a with suitable coefficients. The 
conformal weights of L_4a and L_4(T are K,A = 4 + |. The corre- 
sponding RG eigenvalue is Uo-^^ = — 2 — |. The conformal weight 
of L_^L_2,(^ is /ia.,33 = 6 + |. The corresponding RG eigenvalue is 

• Operators belonging to [e] 

The most important contribution from the [e] family is the one 
proportional to hft which is responsible for the hf term which 
appears in Ut as we discussed in the previous section. Besides 
this one, the lowest operators which appear in the [e] family must 
be of the type L_4e or L_4e. In fact the same mechanism which 
allowed us to eliminate the secondary fields of level 2 in the [a] 
family also works for the [e] family. On top of this in the [e] family 
a new null vector appears at level 3, thus allowing us to eliminate 
also all the fields at this level. Keeping also into account the fact 
that the corresponding Ui{hi) functions must start from hf we 
immediately see that all these operators have too high powers of 
hi to contribute to the scaling function and can be neglected. 

B] Operators which contain L„i,L_i generators. 

The lowest operators are, in order of increasing weight: 
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• L_iL_i(T, whose conformal weight is h^^ii = 2 + |. The corre- 
sponding RG eigenvalue is ya,u = ~|- 

• L^^L^^cr, whose conformal weight is /I0.22 = 4 + |. The corre- 
sponding RG eigenvalue is — —2 — |. 

• L_iZ_ie, whose conformal weight is h^ ii — 3. The corresponding 
RG eigenvalue is y^ ii — —1. 

3.4 Scaling functions 

Using the results of the previous section we are now in the position to write 
the expression for the scaling functions keeping all the corrections up to the 
order hf. 

3.4.1 The free energy 

Due to translational invariance, only the secondary fields which are not gen- 
erated by L_i,L_i contribute to the free energy. We find, for the singular 
part of the lattice free energy: 

fs{hi) = 4 |/ii|Tf(l + A^_J/i^|Tl + A^^2l^d^ +4,3l^i|^ +4,4l^d^ + 

4^5|/i,|i + 4^6l^# ) (82) 

where denotes the amplitude, normalized to the critical amplitude, of 
the n*'^ subleading correction. 

Let us discuss the origin of the various corrections: 

A 1 IT- |i5 

• Af^i|ni|i5 

this term is entirely due to the TT, and irrelevant fields in the 
Hamiltonian. 

1 1 22 

• A{._2|hi|^ 

this term is due to the bthf term in ut (or, equivalently, to the appear- 
ance of a h^e term in the Hamiltonian). 

A 1 111^ 

• Af 3 1 hi 1 15 

this term is due to the ehhf term in Uh. 
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A 1 II \^ 

• AJ- 4 1 hi 1 15 

this term keeps into account the second term in the Taylor expansion 
of the TT like corrections and the contribution of the fields L^^L^^I 
and hL_4^a in the Hamiltonian. 

• A[.g|hi|ft 

this is the product of the Af^i and Af^ corrections. 

• A^6|hi|T5 

this is the second term in the Taylor expansion of the hfe correction. 
To these terms we must then add the bulk contributions 

fb{hi) = U + U,^h^ + U,2ht + --- . (83) 

We have already noticed in sect. 2.1.2 that can be obtained from the 
exact solution of the Ising model on the lattice at the critical point. Also the 
next term: f^^i can be evaluated (with a precision of ten digits) by noticing 
that it corresponds to the constant contribution to the susceptibility at the 
critical point. This term has been evaluated in |2^. We neglect for the 
moment this information and keep the fb,i amplitude in the scaling function 
as a free parameter. It is the first subleading term in the scaling function and 
as such it can be rather precisely estimated with the fitting procedure that 
we shall discuss below. We shall compare our estimates with the expected 
value in sect. 6 and use the comparison as a test of the reliability of our 
results. 

Combining eqs. 

f{hi) = fb + A^f \hi\^'{l + A^f^^\hi\^ + A^jjhil^ + A^j^^lhifl + 

W I , I 30 w I , I 32 w I , I 38 w I , I 44 . /„,s 

A^f^^\hi\-s + A^fJhi\-o + A^fJhi\-o + A^f^Q\hi\-5 ) (84) 

where Ajj, is ^ and A^jq takes also into account now the contribution 

of fb,2- 

Deriving this expression with respect to hi we obtain the scaling functions 
for the magnetization and the susceptibility^. 

"^In this way we obtain directly the lattice definitions of these two quantities, since we 
are deriving the lattice free energy with respect to the lattice magnetic field. There is no 
need to go through the continuum definition of the magnetization. 



82) and (K3) we find: 
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3.4.2 Mass spectrum 

The simplest way to deal with the mass spectrum is to fit the square of 
the masses. The scaling function turns out to be very similar to that which 
describes the singular part of the free energy eq. (P^. The only additional 
terms are due to the secondary fields which contain It turns out 

that the corresponding scaling dimension exactly match those which already 
appear in eq. (^21). In fact 

• /i/L_iL_icr gives a contribution which scales with \hi\^ and its ampli- 
tude can be absorbed in v4j i. 

• hiL'^^L'^ia gives a contribution which scales with \hi\^ and its ampli- 
tude can be absorbed in ^^^4. 

_ie gives a contribution which scales with |/i;|i5 and its ampli- 
tude can be absorbed in ^4^ 5. 

Thus the functional form of the scaling function for the masses is exactly 
the same of eq. (p2|). 

Aljhi\i^ + AlJhi\i^o ) 

However we shall see below that the presence of these new fields and in 
particular of L_iL_i(T has very important consequences. 

3.4.3 Internal energy 

We may obtain the internal energy as a derivative with respect to t of the 
singular part of the free energy. However in doing this we must resume (as 
discussed above) the t dependence in the scaling variables. This leads to some 
new terms in the scaling function with powers ^ (due to the Cht term in Uh), 
II and II (due to the t terms in scaling variables of the irrelevant operators). 
It is nice to see that the presence of these additional contributions can be 
understood in another, equivalent, way. Looking at eq. ([T6|) or (|7^) we see 
that the internal energy on the lattice contains a term of type hia. The 
powers listed above are exactly those that we obtain keeping into account 
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the additional hia term in the scahng function. Keeping also into account 
the bulk contribution we end up with the following scahng function. We 
have: 

E{hi) = 1^/1^(1 + ^kil^'l^ + ^E,2\hi\^' + ^Efilhifl + A^E^i^g\hifllog\hi\ + 
A^E,3\hi\^ + ^£,41^/1^ + ^E,5\hifl + A^Efilhil^ + ^Ejlhil^ + ) (86) 

where denotes the amplitude of the hf term in the bulk part of the 
internal energy, A'-^ i^g denotes the amplitude of the hflog\hi\ term discussed 
in sect. 3.1 and the bulk constant term has been already taken into account 
in the definition of E{h). The first correction which appears in the internal 
energy (with amplitude A^ i) is the one with the lowest power of hi among 
all the subleading terms of the various scaling functions this. Its effect on 
the scaling behaviour of the internal energy is very important and it is easily 
observable also in standard Monte Carlo simulations ||13|| . 

3.4.4 Overlaps 

Also in this case we fitted the square of the overlap constants. The scaling 
functions can be obtained with a straightforward application of the argu- 
ments discussed above. Also fields generated by hL_iL_i must be taken into 
account. Moreover, for the overlaps with the internal energy also the ha 
term must be taken into account. We end up with the following result for 
the magnetic overlaps. 

A^F- ,i\hi\^ + ^^F- + ) • (87) 

While for the energy overlaps we have 

+ ) • (88) 
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4 The transfer matrix method 



We computed the mass spectrum and observables by numerical diagonaliza- 
tion of the transfer matrix. The transfer matrix was introduced by Kramers 



and Wannier in 1941. For a discussion of the transfer matrix see e.g. 
refs. [p8t p9| . The starting point is a simple transformation of the Boltzmann 
factor 

exp /3 (^n(^m + hi^an] =T{ui,U2) T{U2,U3) ...T{ulo,Ui) (89) 

\ <n,m> n j 

where = (o"(„q^i), ct(„(,^2), ••• c"(no,Li)) is the spin configuration on the time 
slice no. T is given by 

7'(Mno,^ino+l) = '^(^ino)^'^^ t/(Mno , Mno+l) V {Un^j^xfl'^ (90) 

with 

f^(Mno,Mno+l) = exp ^ Cr(no,ni)Cr(no+l,ni) ^ (91) 

and 

/ Li Li \ 

y[Un^) =exp /3 ^ CT(„(„ni)CT(no,„i+l) + h X! f^{no,ni) ■ (92) 

\ ni=l ni y 

The partition function becomes 

Z = 5] exp(/3 ^ cr„or^ + /i/E^n) = = E^f° (93) 

(T„±l \ <n,m> n J i 

where T is interpreted as a matrix. The time-slice configurations are the 
indices of the matrix. The number of configurations on a time slice is 2^^. 
Therefore the transfer matrix is a 2^^ x 2^^ matrix. By construction the 
transfer matrix is positive and symmetric. The Aj are the eigenvalues of the 
transfer matrix. 



4.1 Computing observables 

Observables that are defined on a single time slice can be easily expressed in 
the transfer matrix formalism. Let us discuss as examples the magnetisation 
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(94) 



and the internal energy. 

<a,,> = ^ 1 

tr ^T^o _ < i\S\i > 

tr T^o ~ Af ° 

where S" is a diagonal matrix. The values on the diagonal are given by ai^i 
on the configurations. {S{u,u') = 5{u,u') u^^^ where u^^^ denotes a on the 
first site of the time slice). The \i > are normalized eigenvectors of T. 
In the limit Lq ^ oo the expression simplifies to 

< >=< 0\S\0 > (95) 

where |0 > is the eigenvector with the largest eigenvalue. 

The energy can be computed in a similar way. The diagonal matrix 
corresponding to the energy is given by 

E{u,u') =6{u,u') u^^^u^^^ . (96) 

Note that we can only express the product of nearest neighbour spins in this 
simple form if both spins belong to the same time slice. 

In order to understand the relation of the mass spectrum with the eigen- 
value spectrum of the transfer matrix we have to compute correlation func- 
tions with separation in time direction. The time-slice correlation function 
eq. ( P^ ) becomes in the limit Lq — >■ cxo 

< So Sr >= J2 exp(-m, |r|) < 0|^|i > < i\S\0 > (97) 




mi = -logM (98) 



with 



and S = 6{u,u') J^m Note that S is translational invariant (in the 
space direction) and has therefore only overlaps with zero-moment eigenvec- 
tors of T. 

With eq. (|23|) we get 



< 0\S\i > 



rrii Li U . (99) 

< 01510 > ^ ^ 



An analogous result can be obtained for the energy. 
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4.2 Computing the eigenvectors and eigenvalues of T 

The remaining problem is to compute (numerically) eigenvectors and eigen- 
values of the transfer matrix. Since we are interested in the thermodynamic 
limit as well as in the continuum limit we would like to use as large values of 
Li as possible. This soon becomes a very difficult task since the dimension 
of the transfer matrix increases exponentially with Li. The problem slightly 
simplifies if one is interested in the computation of the the leading eigenval- 
ues and eigenvectors only, and in these last years various methods have been 
developed to address this task (for a comprehensive discussion of existing 
approaches see e.g. ref. ||2^ or the appendix of ref. |^). In particular there 



are two approaches which have shown to be the most effective ones. 

• The first one reduces the numerical complexity of the problem by writ- 
ing the transfer matrix as a product of sparse matrices. See refs. 
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• The second one is to reduce the dimension of the transfer matrix by 
restricting it to definite channels. 

Since we are only interested in the zero-momentum states of the system we 
decided to follow the second approach and to compute the zero-momentum 
reduced transfer matrix. The zero-momentum reduced transfer matrix acts 
on the space of equivalence classes of configurations on slices that transform 
into each other by translations. 

The matrix elements of the reduced transfer-matrix are given by 

f{u,v) = {n{u) n{d))-'/^ 5]^T(M,t;) = {n{u)/n{v)y/^Y.T{u,v) . 

(100) 

where n{u) is the number of configurations in u. For example for Lq = 20 
the dimension of the transfer matrix is reduced from 1048576 to 52488. 

Still the matrix is too large to save all elements of the matrix in the mem- 
ory of the computer. Therefore we applied an iterative solver and computed 
the elements of T whenever they were needed. 

As solver we used a generalized power method as discussed in the ap- 



pendix of ref. pO| 
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The lattice sizes that we could reach in this way were large enough for 
our purpose, thus we made no further effort to improve our method and it is 
well possible that our algorithm might still not be the optimal one. 

We propose here, as a suggestion to the interested reader, some directions 
in which it could be improved. 

• One could try to mix the two strategies mentioned above and try to 
factorize the reduced transfer matrix as a product of sparse matrices. 
However note that the complexity of the problem increases exponen- 
tially with the lattice size. Therefore even a big improvement in the 
method would allow just to go up in the maximal Li by a few sites. 

• One could study the transfer matrix along the diagonals of the square 
lattice. Since the distance between two points on the diagonal is ^/2, 
naively one could increase the accessible lattice size by a factor of \/2. 



5 Thermodynamic limit 

In order to take the thermodynamic limit we must know the finite size scaling 
behaviour of the various observables as a function of Li. This is a very inter- 
esting subject in itself and several exact results have been obtained in this 
context starting from the exact S-matrix solution and using Thermodynamic 
Bethe Ansatz (TBA) techniques |31[ . 



For instance, it is possible to construct a large L asymptotic expansion for 
the finite size scaling (FSS) of the energy levels based only on the knowledge 
of the exact S-matrix of the theory |^2|, Let us look to this FSS behaviour 



in more detail. 

Let us define Ama(L) as the deviation of the mass rria of the particle a 
from its asymptotic value: 

AniaiL) = ma{L) - niaioo) . (101) 

Then in the large L limit, the shift (normalized to the lowest mass mi) is 
dominated by an exponential decrease of the type 

^ ~ E ^ exp i-fiabcL) 102 

nil 8nii Habc 
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where the constants jj^abc and Xabc can be obtained from the S-matrix and 
the prime in the sum of eq. (|102|) means that the sum must be done only on 
those combinations of indices that fulfill the condition: — < m^. In 
particular the ^abc turn out to be of order one, so that the FSS corrections are 
dominated (as one could naively expect) by a decreasing exponential of the 
type exp(— Li/^) where the correlation length ^ is the inverse of the lowest 
mass of the theory. 

In principle we could use our data to test also the TBA predictions for the 
FSS. However we preferred to follow a different approach. We chose values of 
h large enough so as to fulfill the condition Li/^ » 1 for the largest values 
of Li that we could reach. In this way we could essentially neglect all the 
details of the FSS functions and approximate them with a single exponential 
(or, in some cases, with a pair of exponentials). In order to study the FSS 
functions one should choose smaller values of h. We plan to address this 
issue in a forthcoming paper. With our choice of h we drastically simplify 
the FSS problem, however nothing is obtained for free. The price we have to 
pay following this route is that we need to know several terms in the scaling 
functions to fit such large values of h. This explains the major effort that we 
devoted to this issue in sect. 3. 

5.1 Numerical extrapolation 

According to the above discussion, for the extrapolation of our data to the 
thermodynamic limit we made no use of the quantitative theoretical results. 
We made only use of the qualitative result that the corrections due to the 
finite Li vanish exponentially. 

We used as ansatz for the extrapolation either 

A{Li) = A{oo) + ci exp{-Li/zi) (103) 

or 

A{Li) = A{oo) + ci exp(-Li/zi) + C2 exp(-Li/2;2) (104) 

where A represents any of the quantities that we have studied. 

In order to compute the free parameters A{oo), ci and zi or A{oo), ci, zi, 
C2 and Z2 we solved numerically the system of equations that results from the 

lattice sizes Ll^maxj -^l,m,ax 1 and Ll^max 2 or Ll^maxj Li^max 1) Li jYiax 2, 
Ll^max ~ 3 and Li max ~ 4. 
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The error of A{oo) was estimated by comparing results where Li^max is 
the largest lattice size that is available and from L[ = Li^max — 1- Mostly 
ansatz ( p.U4| ) was used to obtain the final result. In some cases however the 
numerical accuracy was not sufficient to resolve the second exponential term. 
Then the final result was taken from the ansatz ( |103| ). 



T 


/ 


/, eq. (^) 


/, eq. (^) 


4 





993343441146 














5 





9923^ 


M038449 














D 





992160642059 





992092835647 








7 





992102865951 





9920^ 


^2710940 








8 





9920^ 


^6845804 





9920^ 


^0699493 





9920^ 


^0279141 


9 





9920^ 


^2188258 





9920^ 


^0279123 





9920^ 


^0180502 


10 





9920^ 


^0787928 





9920^ 


^0185903 





9920^ 


^0161487 


11 





9920^ 


^0356320 





9920^ 


^0164020 





9920^ 


^0157709 


12 





9920^ 


^0220728 





9920^ 


^0158619 





9920^ 


^0156931 


13 





9920^ 


^0177480 





9920^ 


^0157225 





9920^ 


^0156758 


14 





9920^ 


^0163514 





9920^ 


^0156853 





9920^ 


^0156721 


15 





9920^ 


^0158958 





9920^ 


^0156752 





9920^ 


^0156716 


16 





9920^ 


^0157458 





9920^ 


^0156722 





9920^ 


^0156709 


17 





9920^ 


^0156961 





9920^ 


^0156715 





9920^ 


^0156713 


18 





9920^ 


^0156795 





9920^ 


^0156712 





9920^ 


^0156710 


19 





9920^ 


^0156739 





9920^ 


^0156710 





9920^ 


^0156710 


20 





9920^ 


^0156721 





9920^ 


^0156712 





9920^ 


^0156713 


21 





9920^ 


^0156714 





9920^ 


^0156710 









Table 6: Extrapolation of the free energy at hi = 0.075 to the thermodynamic 
hmit. In the first column we give the lattice size Li . In the second column 
the free energy for this lattice size is given. In the third column we present 
the extrapolation with a single exponential and in the fourth column the 
extrapolation with a double exponential ansatz. 

In tab. ^ we give, as example, the extrapolation to the thermodynamic 
limit of the free energy at hi = 0.075. As input for the extrapolation we used 
the free energy computed up to 12 digits. We consider all these digits save 
of rounding errors. Within the given precision the free energy has not yet 
converged at Li = 21. The single exponential extrapolation ( |103| ) converges 
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(within the given precision) at Li = 18. For larger lattices the result fluctu- 
ates in the last digit due to rounding errors of the input data. The double 
exponential extrapolation (JTUBj) converges at Li = 16. As final result for the 
thermodynamic limit we quote /(0.075) = 0.99208015671. 



6 Analysis of the results 

The major problem in extracting the continuum limit results from the data 
listed in tabs. p!0| - |T3| is to estimate the systematic errors involved in the trun- 
cation of the scaling functions that we use in the fits. We shall devote the 
first part of this section to a detailed description of the procedure that we 
followed to estimate this uncertainty. We shall give upper and lower bounds 
for the critical amplitudes which turn out to be very near to each other and 
allow for high precision predictions (in some cases we can fix 5 or even 6 sig- 
nificative digits). We then compare our predictions with the results obtained 
in the framework of the S-matrix approach. In all cases we find a perfect 
agreement within our bounds. Finally in sect. 6.3, we give, assuming as fixed 
input the S-matrix predictions for the critical amplitudes, our best estimates 
for the amplitudes of some of the subleading terms involved in the fits. 



6.1 Systematic errors 

In order to estimate the systematic errors involved in our estimates of the 
critical amplitudes we performed for each observable several independent 
fits starting with a fitting function containing only the dominant scaling 
dimension and then adding the subleading fields one by one. For each fitting 



function we tried first to fit all the exiting data (those listed in tabs. p!0|-p 
and then eliminated the data one by one starting from the farthest from the 
critical point (i.e. from those with the highest values of hi). Among the 
(very large) set of estimates of the critical amplitudes we selected only those 
fulfilling the following requirements: 

1] The reduced of the fit must be of order unity f\. In order to fix 
precisely a threshold we required the fit to have a confidence level larger 



^This is a slightly incorrect use of the function since the input data are affected by 
errors which are of systematic more than statistic nature. Notice however that we do not 
use it to determine best fit values for the observables that we fit (we shall only give upper 
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than 30%. 



2] The number of degrees of freedom of the fit (i.e. the number of data 
fitted minus the number of free parameters in the fitting function) must 
be larger than 3. 

3] For all the subleading fields included in the fitting function, the am- 
plitude estimated from the fit must be larger than the corresponding 
errors, otherwise the field is eliminated from the fit. 

4] The amplitudes of the subleading fields (in units of the critical ampli- 
tude) must be such that when multiplied for the corresponding power 
of hi, (for the largest value of hi involved in the fit) must give a contri- 
bution much smaller than 1 (in order to fix a threshold we required it 
to be strictly smaller than 0.3). 

In general only a small number of combinations of data and degrees of free- 
dom fulfills simultaneously all these requirements. Among all the correspond- 
ing estimates of the critical amplitude we then select the smallest and the 
largest ones as lower and upper bounds^. 



6.2 Critical amplitudes 

In tab. 1^ we report as an example the fits to the magnetization (with the 
scaling function obtained by deriving eq. (P^)) fulfilling the above require- 
ments. For each value of Nf we only report the fits with the minimum and 
maximum allowed number of d.o.f., since the best fit result for changes 
monotonically as the data are eliminated from the fit. This is a general pat- 
tern for all the observables that we studied and greatly simplifies the analysis 
of the data. Looking at the table one can see that at least four parameters 
are needed in the fit to have a reasonable confidence level, due to the very 
small error of the data that we use. In the last line we report the only fit in 



which all the 25 data reported in tab. |10| have been used. It required taking 



into account the first eight terms of the scaling function. For Nf > 8, even 

and lower bounds for them) but only as a tool to eliminate those situations in which the 
fitting functions are clearly unable to describe the input data. 

^In making this choice we also keep into account the errors in the estimates induced 
by the systematic errors of the input data. 
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if we use all the data at our disposal we cannot fulfil requirement 3. It is 
interesting to notice that the fits which give the best approximations to the 
exact value of A^j^j are those in which we use the largest possible number of 
terms of the scaling function. This is a general pattern for all the observables 
that we studied. All the fits were performed using the double precision NAG 
routine G02DAF. The bounds that we obtained are listed in tab. ^ together 
with the S-matrix predictions. From these results we immediately obtain the 
upper and lower bounds for the universal amplitude ratios of tab. |[ They 
are reported in the last two lines of tab. |^. 





Nf 


d.o.f. 


C.L. 


1.05898893(196) 


4 


4 


83% 


1.05899447(58) 


4 


6 


50% 


1.05898584(74) 


5 


6 


94% 


1.05898882(22) 


5 


7 


48% 


1.05898178(156) 


6 


6 


98% 


1.05898375(8) 


6 


10 


98% 


1.05898433(8) 


7 


14 


80% 


1.05898694(18) 


8 


15 


99% 


1.05898729(13) 


8 


17 


96% 



Table 7: Fits to the magnetization fulfilling the requirements 1-4 (see text). 
In the first column the best fit results for the critical amplitude (with in 
parenthesis the error induced by the systematic errors of the input data), in 
the second column the number of parameters in the fit, in the third column 
the number of degrees of freedom and in the last column the confidence level. 
For each value of Nf we only report the fits with the minimum and maximum 
allowed number of d.o.f, since the best fit of A^j changes monotonically as 
the data are eliminated from the fit. 



6.3 Subleading operators 

In principle we could try to estimate in the fits discussed above also the 
amplitudes of the first two or three subleading terms in the scaling functions, 
however it is clear that the results that we would obtain would be strongly 
cross correlated and we would not be able to give reliable estimates for the 
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(Jbservable 


Lower bound 


Upper bound 


1 neory 




U.yyz / yoD 


u.yyzouuD 


n 000700C; 
U.yyz /yyo... 


A^ 


1 058980 


1 058995 


1 058986 




n 07055 


07072 


0705Q9 




0.58050 


0.58059 


0.58051... 




4.01031 


4.01052 


4.01040... 


Al 

1712 


6.486 


6.491 


6.4890... 


Al 

1713 


7.91 


8.02 


7.9769... 


I^Ff 1 


0.6405 


0.6411 


0.6409... 


14' 1 


3.699 


3.714 


3.7066... 


/ 1 


0.3 


0.35 


0.3387... 


14' 1 


3.32 


~ 3.45 


3.4222... 




6.7774 


6.7789 


6.77828... 


^2 


3.2296 


3.2374 


3.23514... 



Table 8: howev and upper bounds for various critical amplitudes discussed 
in the text and, in the last two lines, for the two universal amplitude ratios 
and Q2- 

corresponding errors (except, at most, for the first one of them, the next to 
leading term in the scaling function). 

In order to obtain some information on the subleading terms we decided 
to follow another route. The results of the previous section strongly support 
the correctness of the S-matrix predictions. We decided then to assume 
these predictions as an input of our analysis, fixing their values in the scaling 
functions. Then we used the same procedure discussed in sect. 6.1 to identify 
the amplitude of the first subleading field. Let us look to the various scaling 
functions in more detail 

6.3.1 Free energy 

This is the case for which we have the most precise data. Moreover we may 
use the data for the magnetization and the susceptibility as a cross check of 
our estimates. 

Combining all the data at our disposal we end up with a rather precise 
estimate for Ai ^, which turns out to be bounded by: 
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-0.055 < < - 0.050 . (105) 

As mentioned in sect. 3.4.1 it is possible to evaluate this amplitude in a 
completely different way, by looking at the constant term in the magnetic 
susceptibility of the model at the critical point. The comparison between 
our estimate and the expected value represents a test of the reliability of our 



fitting procedure. The expected value of this amplitude ||26[ is (in our units) 

A^^^ = - 0.0524442... (106) 
which is indeed in perfect agreement with our estimates. 



We can then use the value of eq. ( p.06| ) as a fixed input and try to estimate 
the amplitude of the following subleading field which has a very important 
physical meaning being the contribution due to the presence of the TT (and 
related terms) operator in the lattice Hamiltonian. Remarkably enough, it 
turns out, by applying the usual analysis, that the corresponding amphtude 

I is compatible with zero. More precisely we see that, changing the number 
of input data and of parameters in the scaling function, the sign of ^ 
changes randomly and its modulus is never larger than 10~^. The same 
pattern is reproduced in the magnetization and in the susceptibility. We 
summarize these observations with the following bound 

< 0.00005 . (107) 

This result agrees with the observation concerning the absence of corrections 
due to irrelevant operators in the case t ^ and h = 0. (For a thorough 
discussion of this point see and refs. therein.) 



If we also assume that Ajr^ = then we may give a reliable estimate for 
the amplitude of 2 which turns out to be bounded by: 

0.020 < v4'^2 < 0.022 . (108) 

This is the highest subleading term that we could study with a reasonable 
degree of confidence in our scaling functions. 

6.3.2 Internal energy 

In the case of the internal energy the first subleading amplitude can be stud- 
ied with very high confidence since it is associated to a very small exponent: 
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|/i/|i5. The result turns out to be 

- 0.646 < A^E,i < - 0.644 . (109) 

In this case the fits are so constrained that we can study with a rather good 
degree of confidence also the next subleading correction, 2 which is very 
interesting, since it again contains the TT term discussed above. In agree- 
ment with the previous observations also in this case the amplitude turns 
out to be compatible with zero. More precisely its sign changes randomly as 
the input data are changed in the fits and its modulus can be bounded by: 

|A^2l< 0-005 . (110) 



which is not as strong as the bound of eq. ( |107| ) but clearly goes in the same 
direction. 

6.3.3 Masses 

The most interesting feature of the scaling functions for the masses is that 
there is no analytic term and the first subleading contribution Al^, ^ is the 
exact analogous of the ^ term for the free energy. 

In this case we find a non zero contribution for A^, ^. In particular we 
find the following bounds for the three masses that we studied: 

-0.21 < Al^^, < -0.20 (111) 

-0.48 < A^^^^^ < -0.41 (112) 
-0.65 < Al^.^ -^ < -0.50 . (113) 

In the case of the masses a preferred direction is singled out. Therefore, 
one has to expect that there is a finite overlap with the irrelevant operator 
that breaks the rotational symmetry. Notice that a similar contribution has 
been observed also in the case of the thermal perturbation of the Ising model 
where the authors studied the breaking of rotational invariance in the 



m 



two point correlator (see sect.IV-G of for a discussion of this point). 



Our results on the amplitude of the subleading corrections are summa- 
rized in tab. Bl 
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-0.055 



0.020 
0.646 



< 



< 



< 




< 



< 



< 



< 



< 



- 0.050 
0.00005 
0.022 

- 0.644 
0.005 



-0.21 


< 


^mi,l 


< 


-0.20 


-0.48 


< 


^m2,l 


< 


- 0.41 


-0.65 


< 




< 


-0.50 



Table 9: Lower and upper bounds for the amplitudes of some of the sublead- 
ing corrections. 

7 Conclusions 

The major goal of this paper was to test the S-matrix description proposed 
by Zamolodchikov in [Q] for the 2d Ising model perturbed by a magnetic field. 
To this end we developed some tools and obtained some results which are 
rather interesting in themselves. In particular 

• We improved the standard transfer matrix calculations by implement- 
ing a zero momentum projection which allowed us to drastically reduce 
the dimension of the matrix. 

• We discussed in detail the relationship between continuum and lattice 
observables. 

• By using CFT results at the critical point we constructed the first 7-8 
terms of the scaling functions for various quantities on the lattice. 

We could obtain in this way very precise numerical estimates for several 
critical amplitudes (in some cases with 5 or even 6 significative digits) and 
in all cases we found a perfect agreement between S-matrix predictions and 
lattice results. 

By assuming the S-matrix predictions as an input of our analysis we could 
estimate some of the subleading amplitudes in the scaling functions. In one 
case the value of the subleading amplitude was already known and again we 
found a complete agreement between theoretical prediction and numerical 
estimate. For the remaining ones there is up to our knowledge no theoretical 
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prediction. They are collected in tab. |^ and represent the most interesting 
outcome of our analysis. We leave them as a challenge for theorists working 
in the field. 

Among the others, the most surprising result concerns the TT term which 
turns out to have a negligible amplitude in the scaling functions of the trans- 
lationally invariant observables. It would be nice to understand which is the 
reason of such behaviour. 

Let us conclude by stressing that the techniques that we have developed 
can be easily extended to the case in which a combinations of both thermal 
and magnetic perturbations is present. In this case the exact integrability 
is lost and our numerical methods could help to test new approaches and 
suggest new ideas. 
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Tables of data 



Table 10: Data used in the £ts 



hi 


/ 


M 


E 


0.20 


1.106272538601(1) 


0.934113075978(1) 


0.182495416253(1) 


0.19 


1.096943627061(1) 


0.931644255995(1) 


0.179101587939(1) 


0.18 


1.087640179593(1) 


0.929017517063(1) 


0.175544472125(1) 


0.17 


1.078363862266(1) 


0.926215008782(1) 


0.171809915290(1) 


0.16 


1.069116534998(1) 


0.923215694344(1) 


0.167881687799(1) 


0.15 


1.059900287285(1) 


0.919994540350(1) 


0.163741028380(1) 


0.14 


1.050717483321(1) 


0.916521430645(1) 


0.159366050850(1) 


0.13 


1.041570819851(1) 


0.91275968274(1) 


0.154730958303(1) 


0.12 


1.032463401585(1) 


0.90866397795(1) 


0.149804982192(1) 


0.11 


1.023398841451(1) 


0.90417740232(1) 


0.14455091814(1) 


0.10 


1.014381396853(1) 


0.89922709483(1) 


0.13892305302(1) 


0.09 


1.00541615982(1) 


0.89371763122(1) 


0.13286414108(1) 


0.08 


0.99650933082(1) 


0.88752055778(1) 


0.12630083230(1) 


0.075 


0.99208015671(1) 


0.88411094491(1) 


0.1228010112(1) 


0.066103019026467 


0.98424336850(1) 


0.87741739906(1) 


0.1161548337(1) 


0.055085849188723 


0.97462849835(1) 


0.86771621938(2) 


0.10703505648(2) 


0.05 


0.97022834(1) 


0.86255168(1) 


0.10241966(1) 


0.044068679350978 


0.96513182856(1) 


0.8558157835(1) 


0.096641767(1) 


0.033051509513233 


0.95578360408(2) 


0.840485633(1) 


0.084469355(1) 


0.03 


0.95322656(1) 


0.83533709(5) 


0.0806726(1) 


0.022034339675489 


0.9466343376(2) 


0.81901353(2) 


0.0695436(1) 


0.02 


0.94497330(2) 


0.8139196(1) 


0.0663409(2) 


0.015 


0.9409395(1) 


0.7988985(1) 


0.057595(1) 


0.01 


0.936994(1) 


0.77805(5) 


0.047045(3) 


0.0088137358702 


0.93607461(2) 


0.771605(1) 


0.044149(2) 
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Table 11: Data used in the fits 



1, 

III 


L/ nil 


i//"2 


"1 / . > , 
1 / ///.--J 


0.20 


0.59778522553(1) 


0.37795775263(1) 


0.310888(1) 


0.19 


0. 01388448719(1) 


0.38765653507(1) 


0.318578(1) 


U.io 


U.d3134d/U4/ I [I) 


U.39oio99oo29(l j 


U.32D94U(i ) 


U.i / 


u.DoUo / ozo lyjuyi) 




U.OODU 1 1 [^) 


0.16 


0.67120940172(1) 


0.42228634593(5) 


0.346115(3) 


0.15 


0.69415734924(1) 


0.43618773124(1) 


0.357209(3) 


0.14 


0.71959442645(1) 


0.45161985381(4) 


0.369548(4) 


0.13 


0.74799884641(1) 


0.4688779288(2) 


0.38338(1) 


0.12 


0.77998715416(1) 


0.488342470(1) 


0.3990(1) 


0.11 


0.81637015277(1) 


0.510513817(1) 


0.4168(1) 


0.10 


0.85823913569(5) 


0.5360654(1) 


0.4374(5) 


0.09 


0.9071039295(1) 


0.5659287(6) 


0.4624(5) 


0.08 


0.965123997(1) 


0.60144(1) 


0.492(1) 


0.075 


0.998514180(1) 


0.62189(1) 


0.508(1) 


0.066103019026467 


1.067300500(2) 


0.66405(5) 


0.543(1) 


0.055085849188723 


1.17524158(3) 


0.7305(1) 




0.05 


1.237044(1) 


0.768(1) 




0.044068679350978 


1.322589(6) 


0.82(1) 




0.033051509513233 


1.54057(2) 






0.03 


1.6218(2) 






0.022034339675489 


1.91(1) 







50 



Table 12: Data used in the fits 



hi 


1 rcr 12 
F 1 


1 I7'crl2 
F2 1 


n on 
U.zU 


n oon/1 1 0Qfi7i T ('\\ 
u.zyu4iyoo ( 


U.UooUUyool^i ) 


0.19 


0.29570405694(1) 


0.04039999(1) 


0.18 


0.30107729858(1) 


0.04291078(1) 


0.17 


0.30653975241(1) 


0.04554676(1) 


0.16 


0.31209194307(1) 


0.04831337(1) 


0.15 


0.31773424601(1) 


0.05121641(1) 


0.14 


0.32346684419(1) 


0.05426214(1) 


0.13 


0.3292896717(1) 


0.05745711(3) 


0.12 


0.3352023388(1) 


0.0608082(1) 


0.11 


0.3412040323(4) 


0.0643227(2) 


0.10 


0.3472933781(4) 


0.068008(1) 


0.09 


0.3534682486(5) 


0.07187(1) 


0.08 


0.359725487(1) 


0.0759(1) 


0.075 


0.362883627(1) 


0.0780(2) 


0.066103019026467 


0.368548928(2) 


0.0818(5) 


0.055085849188723 


0.3756378(4) 




0.05 


0.378934(5) 




0.044068679350978 


0.38280(2) 




0.033051509513233 


0.3899(1) 
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Table 13: Data used in the fits 



7 

hi 


1 7~'el2 

\r{\ 


1 7~'el2 


n on 
U.zU 


il. 1 UUUi 1404 1 [i) 


o.U40oU0 ( [0) 


0.19 


11.8448924368(1) 


8.3246112(5) 


0.18 


11.9888157747(1) 


8.6017424(5) 


0.17 


12.1315853880(1) 


8.8774943(5) 


0.16 


12.2729879270(1) 


9.1511600(5) 


0.15 


12.4127893980(1) 


9.422023(1) 


0.14 


12.5507307560(1) 


9.689348(2) 


0.13 


12.6865220830(5) 


9.952360(5) 


0.12 


12.819834783(1) 


10.21022(1) 


0.11 


12.950290902(2) 


10.46202(3) 


0.10 


13.077448185(2) 


10.7067(5) 


0.09 


13.200778543(2) 


10.943(3) 


0.08 


13.31963596(3) 


11.17(1) 


0.075 


13.3771415(1) 


11.28(1) 


0.066103019026467 


13.475815(5) 


11.46(2) 


0.055085849188723 


13.59037(3) 


11.6(5) 


0.05 


13.6398(5) 




0.044068679350978 


13.695(1) 




0.033051509513233 


13.78(1) 





52 



